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ABSTRACT 

We use numerical simulations of a cosmological volume to study the X-ray ionisation and 
04 \ heating of the intergalactic medium by an early population of accreting black holes. By con- 

sidering theoretical limits on the accretion rate and observational constraints from the X-ray 

■ background and faint X-ray source counts, we find that the maximum value of the optical 

depth to Thompson scattering that can be produced using these models is r c ~ 0.17, in agree- 
ment with previous semianalytic results. The redshifted soft X-ray background produced by 
these early sources is important in producing a fully ionised atomic hydrogen in the low den- 
sity intergalactic medium before stellar reionisation at redshift z ~ 6 — 7. As a result stellar 

■ reionisation is characterised by an almost instantaneous "overlap phase" of H II regions. The 
■^-j- \ background also produces a second He II reionisation at about redshift three and maintains 
f^S . the temperature of the intergalactic medium at about 10,000 K even at low redshifts. 

If the spectral energy distribution of these sources has a non-negligible high energy 
power-law component, the luminosity in the soft X-ray band of the "typical" galaxies hosting 
intermediate-mass accreting black holes is maximum at z ~ 15 and is about one or two or- 
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ders of magnitude below the sensitivity limit of the Chandra deep field. We find that about a 
thousand of these sources may be present per square arcmin of the sky, producing potentially 
detectable fluctuations. We also estimate that a few rare objects, not present in our small sim- 
ulated volume, could be luminous enough to be visible in the Chandra deep field. XEUS and 
Constellation-X satellites will be able to detect more of these sources that, if radio loud, could 
be used to study the 21 cm forest in absorption. 
rN \ A signature of an early X-ray preionisation is the production of secondary CMB 

^ ■ anisotropies on small angular scales (< 1 arcmin). We find that in these models the power 

spectrum of temperature fluctuations increases with decreasing angular scale (AT ~ 16/iK at 
~ 1 arcsec scales), while for stellar reionisation scenarios the power decreases on smaller 
scales. We also show that the redshifted 21 cm radiation from neutral hydrogen can be 
marginally detected in emission at redshifts 7 < z < 12. At a redshift of about z ~ 30 a 
stronger and narrower (in redshift space) signal in absorption against the CMB, that is pecu- 
liar to these models, could be detectable. 

Key words: cosmology: theory - methods: numerical 



1 INTRODUCTION 

The WMAP satellite measured an optical de pth to Tho mpson 
scattering of the IGM r e ~ 0.17 ± 0.04 <Kogut et all 120031) . 
The result, implies a much earlier start of reionisation with re- 
spect to the redshift when reionisation was completed (at z rc i ~ 
6, as estimated from the absorption spectra of high redshift 
quasars). This has also been interpreted as one of the stronger 
pieces of evidence for the importance of zero-metallicity (Popula- 
tion III) stars in the early universe te.gJCerl2003d:IW\Tt he & Loeb 
[ 20031: ISomerville ; & Liyidl2003l : ICiardi et alJl2003USokasian etafl 
|2004 IChiu et alj|2003l) and of the formation of the first small 



mass galaxies (with Md m ~ 10 6 — 10 8 Mp), a topic that is 
still under debate (e.g., lHaiman et alJll997l: iRicotti et all 1200 ll 
2002b; Vlachacek et al. 2003). In a previous paper (paper I), 
PJ^oTti^^strike^T2004a^ have argued that, even assuming the 
most favorable properties for Population III, it is difficult or im- 
possible to produce the measured value of r e using UV radiation 
from Population III stars. The main problems with this scenario 
are the negative feedback from ionising radiation, mechanical feed- 
back from SN explosions and the contamination of high density 
regions by detritus from the same stars which produce the ionis- 
ing radiation. Taking into account the metal enrichment from SNe 
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and pair instability SNe the Population III epoch turns out to be 
so short lived that Population III stars will never be able to com- 
plete the reionisation of the IGM. But if a significant fraction of 
Population III stars, instead of exploding as SNe, implode into 
black holes (BHs), then the negative feedback on star formation 
is much reduced and the UV radiation produced by thermonuclear 
reactions in the first stars might be important for reionisation. In 
iRicotti & OstrikeJ J2004b) (paper Ha) we have shown that in such 
a scenario, the X-ray background produced by accretion onto stel- 
lar mass seed black holes would be more effective in producing the 
large r e meas ured by WMAP than the UV radiation from the parent 
stars . iMadau et alJl2004l) have also investigated a scenario in which 
accretion on intermediate-mass black holes produces the large opti- 
cal depth to Thompson scattering measured by WMAP. In their pa- 
per they focus is on the physics of black hole mergers and accretion 
onto seed black holes. Their paper is complementary to our stud- 
ies that, instead, focus on radiative transfer processes. Before the 
WMAP measurement of a large r e , t he ef fect of X-ray preion isation 
has been investigated by Ota i200ll> and Venkat esan et alj 1200 ll) . 
The values of r c that they found were smaller since the H I ionisa- 
tion was produced by secondary photoelectrons that can ionise the 
gas to a maximum of 10% ionisation fraction. We find larger values 
of T e , consistent with WMAP, because of the additional ionisations 
produced by the redshifted soft X-ray background. 

In this paper we complement the semianalytic results pre- 
sented in paper Ha by using hydrodynamic cosmological simula- 
tions that include a recipe for star formation and 3D radiative trans- 
fer for hydrogen and helium ionising radi ation. An a pproximate 
solution of the radiative transfer equations iGnedin & AbefcoOll) 
is used to speed up the calculations and make the coupling of the ra- 
diative transfer and hydrodynamic equations computationally fea- 
sible. In addition, radiative transfer for the optically thin X-ray radi- 
ation and H2 dissociating radiation is solved exactly. Line radiative 
transfer in the H2 Lyman- Werner bands (11.3 < hpv < 13.6 eV) 
is also solved exactly for the volume averaged component of the 
radiation field. Using this method we are able to simulate local and 
global radiative feedback effects of galaxy formation on cosmolog- 
ical scales. 

The results of the cosmological simulations are used to pre- 
dict the distinctive observational signatures of the X-ray preionisa- 
tion scenario compared to stellar reionisation models. We find that, 
if this scenario is correct, a new population of X-ray sources that 
do not have optical counterparts may be detected in the Chandra 
deep field. Interestingly, these sources, that might be high redshift 
AGNs but also galaxies hosting a large number of bright Ultra- 
luminous X-ray sources (ULXs), might have been already observed 
iKoekemoer et all20o3) . 

Simulations of a cosmological volume allow us to construct 
maps for the secondary CMB anisotropies produced after recom- 
bination by the first sources of ionising radiation. We also use the 
simulation outputs and the semianalytic models presented in pa- 
per Ha to calculate the redshifted 21 cm radiation from the IGM 
prior to reionisation in emission or absorption against the CMB. In 
the X-ray preionisation scenario the IGM is only partially ionised at 
redshift z > z Ic i — 6 — 7. Therefore, the 21cm signal from the IGM 
at z < 12 can be detected even though it is rather weak. Note that, 
in stellar reionisation scenarios that can produce the optical depth 
to Thompson scattering observed by WMAP, the expected 21cm 
signal is instead undetectable. Finally we estimate the additional 
high energy background due to the postulated population of high 
redshift X-ray sources and its subsequent signatures [(10 ± 5)% in 
the 2-50 keV bands]. 



This paper is organised as follows: in § |2| we show the re- 
sults of cosmological simulations of X-ray preionisation by mass 
accretion on seed BHs. In § [3] we estimate the number of de- 
tectable X-ray point sources at high redshift. In §[4] we calculate 
the redshifted 21cm signal for one of our simulations in absorp- 
tion/emission against the CMB. In §[5] we compute the amplitude 
of the power spectrum of CMB secondary anisotropies on scales of 
a few arcmin produced by X-ray preionisation and stellar reionisa- 
tion scenarios. We summarise the results in §|6|and we discuss the 
observational signatures of the X-ray preionisation scenario when 
compared to stellar reionisation scenarios. 



2 COSMOLOGICAL SIMULATIONS 

In this section we show the results of cosmological simulations in- 
cluding radiative feedback effects and feedback from SN explo- 
sions. The code has been implemented and used extensively to 
study the formation of the first galaxies. In this study we have run 
four new simulations to study the effects of SN explosions and an 
early X-ray background. The simulations were run on COSMOS, 
an SGI Origin 38000 in DAMPT at Cambridge University. 

We adopt a concordance ACDM cosmological model with pa- 
rameters consi s tent w ith the analyses of Sperge fet alJ l2003h and 
iTegmark et all <2004 : fi m = 0.3, fi A = 0.7, h = 0.7 and 
fib = 0.04. The initial spectrum of perturbations has erg = 0.91 
and ris = 1. The box size is L\, ox = 1 h _1 Mpc (comoving) and 
the grid has N'i ox = 128 3 cells. We achieve a maximum mass res- 
olution of Mom = 3.94 x 10 4 h~ x M and spatial resolution of 
488 h' 1 pc (comoving). We fully resolve the star formation in ob- 
jects within the mass range 5 x 10 5 M < M DM < 10 9 M . All 
simulations start at z — 100 and end at z ~ 8. After this redshift 
we stop the simulations because the simulated volume ceases to be 
a statistically representative fraction of the universe. 

2.1 The Code 

The simulations were performed with the "Softened Lagrangian 
Hydro dynamics" (SLH-P 3 M) code described in detail in Gnedin 
1 1995). The cosmological simulation evolves collisionless DM par- 
ticles, gas, "star-particles" and the radiation field in four frequency 
bands: optically thin radiation, Hi , He I and Hell ionising radi- 
ation. The radiative transfer is treated self-consistently (i.e., cou- 
pled with the ga s dynamics and star f ormation) using the OTVET 
approximation (Gnedin & Abel 200 ll) . Star particles and BHs are 
formed as per equation Q in each resolution element that sinks 
below the spatial resolution of the code. The code adopts a de- 
formable mesh to achieve higher resolution in the dense filaments 
of the large scale structure. We solve the line radiative transfer in 
the H2 Lyman- Werner bands for the background radiation, we in- 
clude the effect of secondary ionisation of H and He by X-rays, 
heating by Lya scattering, detailed H2 chemistry and cooling, and 
the absorbed stellar energy distribution (SE P) of the sources tha t 
depends on the mean UV escape fraction iRicotti et alj|20023) . 
H ere we also inc lud e the effect o f SN explosions using the method 
in lGnedirll 19981) . In lRicotti et aljl2002al) we discussed extensively 
the details of the code and the physics included in the simulation, 
focusing on simulations of the first galaxies. We also performed nu- 
merical convergence studies that are especially crucial in this case. 
High mass resolution is needed because the objects that we want to 
resolve have masses 10 s M < M dm < 10 8 Mq. Moreover, the 
box size has to be large enough in order to include at least a few 
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Table 1. Input parameters for hydrodynamic simulations with radiative transfer. 



# 


RUN 


AW 


^box 


Mass Res. 


Res. 


9v 


6* 


eijv(/csc) 


e qso 


Zott 


F\MF 








Mpc 


h' 1 M Q 


h^ 1 pc 














1 


M-PIS 


128 


1.0 


3.94 X 10 4 


488 


Pop III 


0.1 


2 x lO" 3 


2 x lO -3 


8 


1 


2 


M-SN1 


128 


1.0 


3.94 x 10 4 


488 


Pop III 


0.1 


3 x 10~ 4 


2 x 10~ 3 


10 


1 


3 


M-SN2 


128 


1.0 


3.94 x 10 4 


488 


Pop III 


0.1 


3 x 10~ 5 


2 x 10~ 3 


11 


0.1 


4 


M-BH 


128 


1.0 


3.94 x 10 4 


488 


Pop III 


0.1 


3 x 10" 5 


see Eg.l2l 


13 


0.1 



Parameter description. Numerical parameters: NS is the number of grid cells, Lb ox is the box size in comoving h 1 Mpc. Physical parameters: g v is the 
normalised SED (Population III). The meaning of the other parameters is explained in the text (§ 12. 11 . 




Figure 1. Slices through the volume of the simulation M-BH (cf., table0 at redshift z = 19, 12.5 and 9 from left to right. The top panels show the hydrogen 
ionisation fraction (in logarithmic scale) and the bottom panels the gas overdensity (in logarithmic scale). Note the ionisation fraction in the voids is larger 
than in the denser filaments. The overdense gas is smoothed on the filtering scale (i.e., time averaged Jeans length) by the X-ray reheating that precedes the 
IGM partial preionisation. 



of the rare first objects. The first small-mass galaxies will form at 
z ~ 30 — 40 from 3a or more rare density perturbations; the first 
normal galaxies with A/d m 5 x 10 8 M form at z Si 15, also 
from 3cr perturbations. Below, we summarise the meaning of the 
free parameters in the simulations. 

• e„: Efficiency of star formation for the adopted star formation 
law, 



where p« and p gas are the stellar and gas density, respectively, t* is 
the maximum between the dynamical and cooling time. 

• euv : Energy in ionising photons emitted by stars per rest mass 
energy of H atoms (muc 2 ) transformed into stars. This parameter 
depends on the IMF and stellar metallicity. We use Population III 
stars SED. 

• (/esc): escape fraction of ionising photons from a cell. It is 
resolution dependent. 



• e qso : Energy in ionising photons emitted by quasars per rest 
mass energy of H atoms (mnc ) transformed into stars. This means 
that, if e qso = const, the total UV luminosity from quasars is pro- 
portional to that from massive stars. We use the template spectrum 
given in § 3.1 of paper Ha. The spectrum is a double power law 
with a soft X-ray bump produced by absorption of UV photons by 
the obscuring torus and the ISM around the BH (in particular the 
absorption is produced by an H I column density of 10 19 cm -2 and 
He with ionisation fractions xnti/xm = 0.1, a;Hcii/a:Hi = 0.1). 
In paper I we mentioned t hat our spectrum is similar to the tem- 
plate spectrum derived by Sazonov et al. 1 2004). This statement is 
actually misleading since our spectrum is softer, having a cutoff at 
about 100 eV instead of 1 keV. We will discuss in t) !2.3l the impor- 
tance of soft X-rays and we will quantify how the cutoff energy in 
the quasar spectrum affects our results. The energy density in the 
X-ray bands is about /3 = 20% of the total and 40 % of the energy 
of H ionising radiation ( ex = 0.4e qso , because of the local ab- 
sorption of UV photons). The accretion rate at the Eddington limit 
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with efficiency e = 0.2, defined in § 3 of paper Ha, is given by 
0.2/3p ac = exp*, where p* is the star formation rate and p ac is the 
accretion rate onto BHs. 

• z ff : redshift at which the X-ray emission by quasars turns off. 
We assume that e qso = at z < z «. 

• Fimf : This parameter is proportional to the mean metallicity 
yield and energy input by SN explosions of the stellar population. 
We have Fimf = 1 for Population II stars with a Salpeter IMF. But 
Population III stars could have Fimf > 1 if the IMF is top-heavy 
and pair-instability SNe are dominant. If the IMF is dominated by 
subluminous BH forming SNe, or by stars with M* > 260 Mo 
that collapse directly into BHs without exploding as SNe (see § 3 
of paper Ha), we have Fimf < 1. 

We have run four simulations that differ mainly in the epoch at 
which the early accretion onto seed BHs takes place (cf, table0. 
In paper Ha we have studied three models: an early preionisation 
model (with preionisation starting at redshift z pro ~ 25), an in- 
termediate preionisation model (z pro ~ 20) and late preionisation 
model (z prc ~ 15). Due to the limited volume and resolution of 
our simulations, the first sources form at redshift z ~ 27, later than 
in the semianalytic models that do not suffer of these limitations. 
The radiation background from the first sources (cf., Fig.|6j builds 
up to relevant values only about one Hubble time after the forma- 
tion of the first source, corresponding to a redshift z ~ 20. At 
this redshift the ionisation fraction of the IGM will start increasing. 
For this reason, using numerical simulations, we cannot simulate 
the early preionisation model presented in paper Ha. In order to 
simulate earlier object formation, we would need simulations with 
comparable mass resolution to the present but of a larger volume 
of the universe. In such a simulation the first galaxies would have 
formed at redshift z ~ 40 from rare (i.e. , 5<r) peaks of the Gaussian 
perturbations in the density field, allowing the earlier production of 
seed BHs and preionisation by the X-ray background. 

In the first three simulations in table[T](M-PIS, M-SN1 and M- 
SN2) we have assumed a step function for the efficiency of X-ray 
emission e qso (i.e., e qso = 0.002 at z > z g and zero afterwards). 
This means that in the first three simulations we assume that the 
local black hole accretion rate is proportional to the local SFR and 
we explore the effects of changing the stellar IMF and (/esc) - In the 
simulation M-PIS the IMF is top-heavy, {/ e8 c) = 1 and the energy 
input from SN explosions is 7 times larger than for a Salpeter IMF. 
The parameters of this first simulation describe a scenario in which 
Population III stars are supermassive and a substantial fraction of 
them explodes as pair-instability supernovae. In the simulation M- 
SN1 we assume a Salpeter IMF and (/esc) = 1. The parameters of 
this simulation are also consistent with a top-heavy IMF combined 
with (/esc) < 1. The simulation M-SN2 is the same as M-SN1 
but has (/esc) = 10% and the effect of SN explosions is 10 times 
smaller. The parameters of this simulation are also consistent with 
a mildly top-heavy IMF combined with (/esc) < 10% and sublu- 
minous or negligible SN explosions. This assumption is justified by 
the fact that zero-metallicity stars, if more massive than ~ 260 Mq, 
may collapse directly into BHs without exploding as SN. Another 
possibility is that the energy of SN explosions in zero-metallicity 
stars with masse s ^ 100 Ma, is smaller than the canonical value 
F = 10 51 ergs <Umeda & Nomotol2003h . 

The fourth simulation (M-BH) is the same as M-SN2 but has 
larger and time-dependent efficiency of X-ray emission, e qso . The 
X-ray emissivity is similar to the one used in the semianalytic 
model M3 (the intermediate preionisation model) of paper Ha and 
it is physically motivated in § 3 of that paper. The time dependent 



function for e qso is given by 



0.2 x 



exp[(39/(l + Z)) 1 - 5 ] if z > 17 
0.0016((1 + z)/15) 20 ifz< 17. 



(2) 



The efficiency e qso has a maximum at 2 = 17 and then quickly 
decreases. This is the most interesting simulation because the func- 
tion e qso is a fit to a realistic model for accretion and because, as 
shown in the next section, produces a value of the optical depth in 
agreement with th WMAP measures (cf., § 12.21 . 



2.2 Results 

The soft X-ray and hard UV flux from the first accreting BHs are 
initially the most efficient in ionising the intergalactic gas in the 
immediate vicinity of each source. Thus the topology of IGM ioni- 
sation is initially characterised by small Stromgren spheres around 
the most luminous sources. But, due to the long mean free path 
of X-ray photons, the background produced by distant sources 
quickly dominates the ionisation rate and the IGM becomes par- 
tially ionised almost uniformly. Due to the negligible recombi- 
nation rate, the voids have a larger fractional ionisation than the 
denser filaments. The evolution of the topology of preionisation is 
shown in the top panels of Fig. Q where we show the neutral frac- 
tion of the IGM in a slice trough the simulation number four in 
table Q at redshift z — 19, 12.5 and 9. The bottom panels show 
the gas overdensity in the same simulation. It is evident from the 
smooth appearance of the filamentary structure that the X-rays, 
heating the IGM and increasing the Jeans mass, have smoothed the 
gas on scales (i.e., the "filtering" scale) larger than the size of the 
dark matter filaments. 

In Fig.|2fleft) we show the BH accretion rate (top panel) and 
the baryon fraction in BHs, ojbh (bottom panel), for the four simu- 
lations whose parameters are listed in tableQ 

In paper Ha we have discussed the physical requirements 
needed to produce such accretion histories. Here we do not discuss 
this matter further but we note that, given the assumed SED of mini- 
quasars, it is possible to determine an upper limit for the global ac- 
cretion rate. This limit is determined by the requirement that the 
early sources of X-rays contribute to less than 10 — 20% of the 
observed X-ray background in the 2 — 10 keV bands. The sources 
that constitute the bulk of the soft X-ray background has now been 
resolved and identified (mainly Seyfert galaxies at z ~ 1 — 2). Re- 
cently De Luca & Molendi 12004) have estimated that that about 
20% of the observed background may be produced by a new pop- 
ulation of faint X-ray sources, currently undetected within the sen- 
sitivity limits of the deepest X-ray surveys. It does not seem that 
the faint high redshift optically discoverable sources ar e abundant 
enough to contribute much to the background iHunt et alj200l) . So 
there is some room for the population of sources that we are postu- 
lating to exist at high redshift. The solid line in Fig.|2jleft) shows 
the value of the global accretion rate as a function of redshift that 
would produce the 20% of the X-ray background at z — that it is 
still unaccounted for by lower redshift AGNs. 

In Fig. fright) we show the star formation rate (top panel) 
and the baryon fraction in stars, lj* (bottom panel), for the models 
in table Q Depending on the simulation parameters, the thermal 
feedback on the IGM or internal feedbacks from galactic winds due 
to photoevaporation and mechanical energy from SN explosions, 
are dominant in reducing t he global star formation rate in small- 
mass galaxies lOstriker & Gnedirll996tlRicotti et all2002r|) . 

The results of the simulations confirm the results found using 
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z+1 z+1 

Figure 2. (left) Black hole accretion history for the four models in tabled (right) Star formation history for the four models in table0 The simulations differ 
in the efficiency of X-ray emission (duty cycle), IMF and ( f CS c). The solid line shows the upper limit of the accretion rate for the SED we have adopted. If the 
global accretion rate on seed BHs is below the solid line, their contribution to the observed X-ray background at z = (at about 10 keV) is less than 20%. 



the semianalytic models presented in paper Ha. The main signa- 
tures of X-ray preionisation are illustrated in Fig. [3] that plot the 
hydrogen ionisation history (top panel) and thermal history of the 
IGM (bottom panel) as a function of redshift for the four simula- 
tions in table Q The early reheating of the IGM to a temperature 
T ~ 10 4 K is reached when the hydrogen ionisation fraction be- 
comes larger than 10 % (cf., Figs. [5}. This happens because the 
ionisation rate from energetic secondary photoelectrons becomes 
inefficient when the fractional ionisation is larger than 10%, con- 
sequently most of the energy of the X-ray photons is deposited 
into heat. If the X-ray background is large enough, the soft X-ray 
photons emitted from distant sources and redshifted into softer UV 
photons, can reionise the voids above the 10% fractional ionisation. 

The IGM optical depth to electron Thompson scattering, r c , 
and the visibility function, 

g(z) =exp(-T e )-^- ; (3) 

where r] is the conformal time, are plotted in Fig. |4]for the simu- 
lations in table Q We find that, for the most extreme model (run 
number four in table 0, X-ray preionisation can produce an op- 
tical depth to Thompson scattering r c = 0.17, in agreement with 
WMAP measurement. The shape of the visibility function deter- 
mines the power spectrum at large angular scales of the polarised 
CMB radiation (EE) and the temperature-polarisation cross corre- 
lation (TE). In particular the EE power spectrum can be used to 
distinguish between models with the same r c but different ionisa- 
tion history because its shape is determined by the redshift at which 
the visibility function has a peak I Holder et al. 2003). As was found 
in paper Ha the intermediate preionisation scenario is the most ef- 
ficient in producing a larger optical depth to Thompson scattering 
because it is the earlier preionisation scenario that requires the min- 
imum number of ionising photons per baryon (i.e., the recombina- 
tions are negligible). 

Contrary to the double reionisation models proposed bv lCenl 



l2003bl) . in the X-ray preionisation models hydrogen is fully 
reionised only once at low redshift (e.g., by Population II stars 
at z ~ 6 — 7). But in most X-ray preionisation models Hell is 
reionised twice. The helium is doubly ionised a first time at high 
redshift. Afterward, in most simulations, due to the fast decline in 
accretion rate onto seed BHs, it partially recombines before red- 
shift nine. But the rate of He III recombination is slow because of 
the photoionisations from the redshifted X-ray background pho- 
tons. This is shown in Fig. [5] where we plot the He ionisation his- 
tory as a function of redshift for the four simulations in table Q In 
most simulations He II becomes about fully ionised (60-80 % frac- 
tional ionisation) before H I reionisation. It is possible to reionise 
He II before fully reionising H I in the voids if it is the background 
that dominates the ionisation rate (note that the background spec- 
trum has almost no UV photons). Indeed the soft UV photons can- 
not ionise the voids because they are absorbed locally while harder 
photons can still ionise H I but less efficiently than He II . 

At redshifts z £ 9, due to the redshifted X-ray background, 
the ionisation fraction of He III starts increasing again and He II be- 
comes raughtly fully ionised a second time at z ~ 3 (see paper Ha 
for a discussion on this effect). The temperature of the IGM, due to 
He II ionisation, remains almost constant from redshift 6 to redshift 
3 at about 10,000 K (cf., paper Ha), in a greement with observation 
of the line widths of the Lyman-a forest iRicotti et all200rj) . 

In Fig.|6|is shown the background radiation at z — 11.5, 7.8 
and the extrapolated value at z = for model M-BH. The observed 
X-ray background radiation at z — is shown with a thick solid 
line. Even though early accretion on seed black holes produces only 
10% of the observed X-ray background (in the 2-10 keV band), at 
redshift z ~ 3 they may have been predominant, indeed most of 
the X-ray background is produced by Seyfert galaxies at z ~ 1. 
Our analysis of the model constraints posed by the observed unre- 
solved component of the X-ray background (presented in paper Ha 
and also shown here in Fig. [6) are consistent with the findings of 
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Figure 3. Hydrogen ionisation history (top panel) and thermal history (bot- 
tom panel) for the four models in table0 In models of reionisation by Pop- 
ulation III stars hydrogen has usually two distinct epoch of reionisation, the 
first at redshift z ~ 17 — 10 and the second at z ~ 6 — 7. In the X-ray 
preionisation models hydrogen is partially reionised at early times but it 
is fully reionised only once at low redshift (i.e., by Population II stars at 
z ~ 6 — 7). Note that the reheating of the gas always precedes the reioni- 
sation. This characteristic is in common with all the reionisation scenarios. 




1+Z 

Figure 5. Helium ionisation history for the four models in tabled In most 
X-ray preionisation models He II is reionised twice. The helium is doubly 
ionised a first time at high redshift. Afterwards it partially recombines be- 
fore redshift 8-9. At redshifts z ^ 8 — 9, due to the redshifted X-ray back- 
ground, the ionisation fraction of He III starts increasing again and He II 
become fully ionised a second time at z ~ 3 (see paper Ha for a discussion 
on this effect). This has also the effect of keeping the temperature of the 
IGM at about 10,000 K from redshift 6 to 3. 
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Figure 4. Thomson scattering optical depth, r e , and visibility function, 
g(z), as a function of redshift for the simulations in tableQ A simulation 
without X-ray preionisation with stellar reionisation by Population II stars 
at z ?s 7 would have r e = 0.6. 



( Diiks tra et alJ2004l) . They also do not rule out a high redshift mini- 
quasar population that could partially reionise the IGM to 50% ion- 
isation fraction. In particular, in the case that the quasar spectrum 
has only a soft X-ray component (e.g., produced by the hot multi- 
colour accretion disk) and no hard X-ray emission, the X-ray back- 
ground does not pose any constraint on the number density and 
luminosity of a putative high redshift population of X-ray sources. 
Sources with this SED could preionise the IGM very efficiently. 

The efficiency of secondary ionisation from energetic photo- 
electrons is large when the ionisation fraction of the IGM is less 
than 10%. This is illustrated in Fig.Qthat shows the mass-weighted 
distribution of the mean hydrogen ionisation, ihu, versus the over- 
density for the M-BH simulation at redshift z ~ 15 (figure on the 
left) and 10 (figure on the right). At redshift z ~ 15 the hydro- 
gen in most of the IGM volume has a fractional ionisation of 10% 
or less as expected. Interestingly, at z ~ 10, even if reionisation 
by stellar sources is not complete, the underdense regions are al- 
ready almost fully ionised. This is due to the intense background in 
the UV bands arising from the redshifted X-rays emitted by distant 
sources that ionises the IGM, preferentially in the low density re- 
gions that occupy most of the volume. The ionisation fraction in the 
voids remains large because the recombination time is longer than 
the Hubble time in these underdense regions. Because the atomic 
hydrogen in the low density intergalactic medium is almost fully 
ionised before stellar reionisation at redshift z ~ 6 — 7, it follows 
that stellar reionisation is characterised by an almost instantaneous 
"overlap phase" of H II regions. 

In the remaining paragraphs of this section we describe the 
results of each simulation in greater detail. 

In the simulation M-PIS (dot dashed line), the IMF is top- 
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overdensity overdensity 

Figure 7. Mass-weighted distribution of the mean hydrogen ionisation fraction versus the overdensity for the M-BH simulation at 2 ~ 15 (left) and z ~ 10 
(right). The stellar sources did not reionise the IGM at this redshift but the underdense regions are already almost fully reionised. UV photons from the 
redshifted X-ray background ionise the IGM, preferentially the low density regions. The voids remain ionised because the recombination time is longer than 
the Hubble time. 
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h p v (eV) 

Figure 6. Background radiation at 2 = 11.5,7.8 and the extrapolated 
value at z = for model M-BH. To guide the eye, the thin vertical lines 
show the energy of the Lyman continuum of H I , He I and He II . The la- 
bels show the Lyman-a lines of H I , He I and He II . The observed X-ray 
background radiation at z = is shown with a thick solid line. Even though 
early accretion on seed black holes produce only 10% of the observed X-ray 
background (in the 2-10 keV band), at redshift z ~ 3 they may have been 
predominant, indeed most of the X-ray background is produced by Seyfert 
galaxies at z ~ 1. 



heavy, (/ esc ) = 1 and the energy input from SN explosions is 7 
times larger than for a Salpeter IMF. This case is one in which 
pair-instability SNe are important in polluting the IGM and in pro- 
ducing strong mechanical feedback in the ISM. The star formation 
is strongly suppressed by the SN explosions. At z ~ 7 the SFR is 
5 x 10~ 3 Mq Mpc _3 s _1 , 20 times smaller than the observed value 
at z ~ 4 — 6. The box size of 1 h _1 Mpc is too small to achieve 
convergence at z ^ 12 (this is demonstrated more quantitatively in 
paper I), therefore at those redshifts the star formation is probably 
underestimated. But the result shows that star formation is strongly 
suppressed in galaxies less massive than 10 9 Mq for this choice 
of the IMF. Metal pollution from SN explosions may suppress the 
formation of Population III stars before z = 10, when we stop our 
simulation. This would reduce the UV emissivity, euv, and the SN 
energy output, Fimf- Here, we did not take into account this ef- 
fect. In this simulation the assumed efficiency of X-ray emission, 
e qS o = 2 x 10 -3 , is too small to have a substantial effect on the 
ionisation and re-heating of the IGM. Reionisation of H I and He II 
happens at z ~ 8 and the gas is heated to T « 3 x 10 4 K. 

In the simulation M-SN1 (dashed line), we assume a Salpeter 
IMF and (/esc) = 1. Here the global SFR is about 10 times larger 
than in the M-PIS simulation, because of the less violent feedback 
from SN explosions. The increased SFR produces a larger X-ray 
emission that at z ~ 12 re-heats the IGM to T ~ 10 4 K and par- 
tially ionises the IGM (x a ) ~ 10%. The X-ray emission continues 
to increase until redshift z = 10 as a result of the increasing star 
formation and the electron fraction reaches (x c ) ~ 50%. At z ~ 9 
He II is reionised. In this model after z ff = 10 we assume zero 
X-ray emissivity (e qso = 0). The He II starts recombining but the 
ionisation fraction continues to increase, but more slowly. This ef- 
fect is due to UV photons arising from redshifted X-rays of the 
background radiation. For this model r c ~ 0.11 and the visibility 
function peaks between redshifts 6 < z < 9. 

The simulation M-SN2 (short dashed line) is the same as M- 
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Table 2. Results of the hydrodynamic simulations with radiative transfer. 



RUN 


To 


-max 

z g 


y-par 




/XRB 


r(Hen)t H 








xlO -7 


xl0~ 5 


(%) 


(%) 


M-PIS 


0.10 


8 


3.0 


0.1 


19(12) 


> 200 


M-SN1 


0.11 


9 


3.1 


0.2 


11(6) 


60 


M-SN2 


0.13 


10.5 


3.3 


0.3 


9(5) 


25 


M-BH 


0.17 


12 


4.4 


2.0 


17(8) 


20 



Meaning of the values in each column: r c is the Thomson scattering optical 
depth; z™ ax is the high-redshift maximum of the visibility function g(z) 
(it also roughly coincides with the redshift of early reionisation of He II 
; y is the Compton distortion parameter; ojbh is mass fraction in BH in 
units of the baryon cosmic density at redshift z = 8 (to translate ubh to 
cosmic BH mass density, Pbh> multiply by 5.5 X 10 9 Mq Mpc - 3); /xrb 
is the fraction of the X-ray background at 50-100 keV (and 2-10 keV, in 
parenthesis) due to early black holes; r(He h)<h ~ (tjj /nu sn )dnue U /dt 
is the fractional rate of He II photoionisation per Hubble time per helium 
atom at redshift z ~ 2 — 3 due to the redshifted X-ray background. 

SN1 but has (/ CS c) = 10% and the effect of SN explosions is 10 
times smaller. This assumption is justified by the fact that zero- 
metallicity stars may collapse directly into BHs without explod- 
ing as SN. Another possibility is that the energy of SN explo- 
sions in zero-metallic ity stars is smaller than the canonical value 
E = 10 51 ergs JUmeda & Nomotol2003h . The results for this sim- 
ulation are analogous to M-SN1 but by virtue of the larger SFR (by 
three times), partial X-ray ionisation and re-heating starts earlier (at 
z ~ 14) and r ~ 0.13. 

Finally, simulation M-BH (solid line) is the same as M-SN2 
but has a larger X-ray emissivity, e qso . Contrary to previous models 
where we have assumed a constant value of e qso , here the emissivity 
is a fit to the semianalytic intermediate preionisation model (model 
M3 in paper Ha). In this simulation the thermal feedback produced 
by the X-ray background on the SFR is evident. At 2 > 15 the 
global SFR is reduced by about one order of magnitude. It starts 
increasing again and reaches the same magnitude as in the M-SN2 
simulation at z = 12 when more massive galaxies start to be nu- 
merous. The X-ray sources are in the most massive galaxies in 
which the star formation rate (and the accretion rate) is not sup- 
pressed by feedback processes. This can be seen also in the bot- 
tom panels of Fig. Q The X-ray background re-heats and partially 
ionises the IGM starting at redshift z ~ 20. At z ~ 13 He II is 
almost completely reionised and (x e ) ~ 0.7. Afterwards, Hell 
slowly recombines but H I remains partially ionised until z — 7, 
when reionisation by stellar sources is completed. 

In table [2] we summarise the values of some relevant quanti- 
ties computed from the simulations in table Q The meaning of the 
quantities shown in each column of the table are given in the foot- 
note. 

2.3 Dependence on the assumed mini-quasars spectra 

We have seen that in our models an important contribution to H I 
and He II reionisation is produced by redshifted soft X-rays. There- 
fore if the quasar spectrum is strongly absorbed we expect that the 
importance of X-rays for the H I preionisation and He II double 
reionisation may be reduced. 

Our template spectrum has a cuto ff at photon en ergies of a 
few 100 eV and therefore differs from Sazonov et al. ( 20040 ab- 
sorbed spectrum that does not have emission below few keV. We 
have assumed that the mini-quasar spectrum is absorbed by a col- 
umn density of H and He with neutral hydrogen column den- 
sity iV a b — 10 19 cm -2 (see § 12. II . substantially smaller than in 
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Figure 8. Dependence of the optical depth to Thompson scattering of the 
IGM, T e , on the average absorbed mini-quasar spectrum. We show t c as 
a function of the column density of neutral hydrogen absorbing the quasar 
spectrum or, equivalently, the cutoff energy of UV and soft X-ray photons 
in the spectrum (shown on the top axis) for model M4 in paper Ha assuming 
Population III star SED (dashed line) or Population II star SED (solid line). 
The results (shown by the points) are calculated using a semianalytic model 
(see paper Ha) that is equivalent to run M-BH in tableQ 

ISazonov et alJ J2004). There is no strong theoretical argument in 
favour of high redshift quasars being strongly absorbed and having 
no soft X-ray emission. On the contrary, high redshift galaxies have 
gas components with radii of a few 100 pc and mean gas densities 
of 10 — 100 cm -3 . The hydrogen column density is therefore only 
N a b ~ (shi/0.001)10 19 cm -3 , where xn, is the hydrogen neutral 
fraction. Most important, the metallicity is expected to be subso- 
lar (Z ~ 0.01 Zq) and in highly absorbed quasa rs most of the 
X-ray absorption is produced by metal lines. The ISazonov et alJ 
(2004) template is therefore more appropriate for quasars at red- 
shifts z £ 5, hosted by quite massive galaxies. 

More quantitatively, our assumption can be justified by the fol- 
lowing arguments. In our simulations the spatial resolution is about 
r = 50 pc therefore, for the assumed column density 7Y a b = nmr, 
we have nm = nxm ~ 0.1 cm -3 , where n is the gas density and 
ihi is the hydrogen neutral fraction. For a gas in photoionisation 
equilibrium we have 

-( M «-'»")(l^)(i)'(wr)' 

(4) 

where So is the number of ionising photons emitted per second 
and 7? is the distance from the source. For example, a 10 4 Mq 
black hole accreting at the Eddington limit would emit about So ~ 
10 50 - 10 51 s" 1 EUV photons. Given our choice of the absorbed 
spectrum we infer a density profile inside our resolution element 
n = (60 cm _3 )(50 pc/i?) (So/10 51 s _1 ) 1/2 , that seems a quite 
reasonable assumption. 

Nevertheless, it is interesting to study how our results would 
change if the spectrum has a cutoff at energies hv CX A > 100 eV. 
We have run a set of semianalytic models assuming that the X-ray 
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sources are absorbed by larger H I column densities ranging from 
N ah = 10 19 cm" 2 to 10 25 cm 2 (approximately we have hu cu t ~ 
100el/(iV a b/10 19 cm" 2 ) 1 ' 3 )- We found that r c decreases linearly 
with increasing log 7V a b and reaches an asymptotic minimum value 
of T e ~ 0.1 for TV > 10 22 cm -2 . In Fig.[8|we show the dependence 
of the optical depth to Thompson scattering of the IGM, r c , on the 
average absorbed mini-quasar spectrum. We can fit the result shown 
by the solid line with the function 

r e = Min[0.177 - 0.026(log /V ab - 19), 0.1]. (5) 

The reionisation of He II at z ~ 13 shown in Fig. [5] is also 
produced by soft X-ray photons emitted by mini-quasars. Increas- 
ing the absorbing column density of the quasar spectrum produces 
a high redshift reionisation only in the low density regions and a 
partial ionisation in the overdense regions of the IGM. When the 
column density is 7V a b > 10 21 — 10 22 cm -2 , the high-redshift 
reionisation does not happen at all. Similarly, the importance of 
redshifted X-rays for the low redshift reionisation decreases and is 
slightly delayed to redshifts z ~ 2 — 3 when the absorbing col- 
umn density increases. Nevertheless, even for large column densi- 
ties (iV ab ~ 10 24 cm" 2 ), redshifted X -rays are still able to reionise 
He II in underdense regions at z ~ 2 and partially ionise the over- 
dense regions. Though, if we assume a Population III SED for the 
stars and H I reionisation at z ~ 6.5, then He II reionisation is dom- 
inated by the stars and happens earlier, at z ~ 3 — 4. Note that we 
do not include the contribution from observed AGNs at redshifts 
z £ 5. Perhaps these AGNs are sufficient to produce He II reioni- 
sation and IGM reheating at z ~ 3 without the need for additional 
contributions. 

Finally, as mentioned earlier, it is plausible that the spectra of 
high redshift X-ray sources have only a soft X-ray component pro- 
duced by thermal emission from a multicolour accretion disk and 
no hard X-ray emission. Accreting intermediate mass black holes, 
because the disk is hotter than in supermassive black holes, emit in 
the FUV and soft X-rays. This would be the most favourable sce- 
nario for achieving large values of r c since soft X-rays efficiently 
pre-ionise the IGM but high redshift sources, even if very numer- 
ous, would be invisible in X-ray and optical deep fields and would 
not contribute to the observed X-ray background. 



3 NUMBER COUNTS OF POINT SOURCES IN THE 
X-RAY BANDS 

We have shown that the high redshift X-ray sources (that we pos- 
tulate to produced the observed r e of the IGM) produce less than 
20% of the observed X-ray background. In this section we calculate 
their contribution to the luminosity function of X-ray point sources. 
Available Chandra deep field observations can already constrain the 
most extreme models (paper Ha, Diikstra et al.| 200^) and the fu ture 
observations with Constellation-X iWhite & Tananbauml 1999b and 
XEUS will be a factor ~ 10 times fainter than the current deep sur- 
vey limit. 

It is important to notice that in these calculations we assume a 
spectral energy distribution of the sources that has a non-negligible 
high energy power-law component (cf., Fig. 4 in paper Ha), in 
agreement with observations of QSOs and ULX in the local uni- 
verse. But if the source spectra are dominated by a multicolour 
disk thermal component their contribution to both the X-ray back- 
ground in the 2-50 keV bands and the faint source counts would 
be negligible. Instead the full ionisation of atomic hydrogen in the 
low density IGM before redshift z ~ 7, the He II reionisation at 
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Figure 9. Number counts of X-ray sources in the 0.5-2 keV band for the 
M-BH run (solid histogram). The brighter tail of the luminosity function 
is produced by the most massive galaxies at redshift z ~ 17 hosting ei- 
ther a supermassive BH of about 10 5 Mq accreting at near the Edding- 
ton rate or several intermediate-mass BHs with the same total mass. The 
dashed portion of the histogram is an extrapolation based on the Press- 
Schechter formalism. The thin vertical line shows the sensitivity limit of 
F x = 2.5 X 10~ 17 erg s _1 cm~ 2 , of the 2 Ms Chandra deep field in the 
0.5-2 keV band, the triangles show the number counts of point sources and 
the solid line a model prediction for the expected counts from AGNs. Fi- 
nally the closed polygon shows the expected steepening of the counts based 
on the statistics of pixel fluctuation from unresolved point sources. The ob- 
served steepening of the number counts slope at fluxes Fx ,$ 10~ 16 is 
consistent with the predictions of our simulation but it is also the predicted 
contribution of starburst galaxies. The high-z sources, contrary to starburst 
galaxies, do not have optical counterparts but could be detectable in the 
infrared if they are not obscured locally by dust. 



z ~ 3 and reheating of the low-redshift IGM is produced by the 
redshifted X-ray background independently of the assumed spectra 
of these early sources. 

In our models we do not make assumptions about the mass 
function of BHs in each galaxy, but only on their total mass and 
accretion efficiency. Two scenarios are therefore plausible for the 
X-ray emission. If seed BHs quickly merge into a massive BH at 
the centre of galaxies through dynamical friction, the X-ray sources 
will look like small AGNs (mini-quasars). If the seed BHs do not 
merge efficiently, in each galaxy we could have several off-centre 
X-ray sources produced by accreting intermediate mass BHs or- 
biting the galaxy. These sources would be similar to ULXs ob- 
served in nearby galaxies, but more numerous (tens or hundreds per 
galaxy) or more luminous. This new type of X-ray sources could 
be common in the early universe if primordial galaxies are com- 
pact and gas rich. This is expected in young galaxies with masses 
Md m Si 10 s M(7) in which stellar winds did not evaporate all their 
gas (Whalen et al. 2003; Ricotti et al. 2004) and star formation did 
not have time to consume most of the gas. Moreover, in the preion- 
isation models, the number density of intermediate mass BHs was 
larger in primordial galaxies than in today's galaxies as was shown 
in Fig. 3 of paper Ha. 



© 2002 RAS, MNRAS OOO.ITWbI 



10 M. Ricotti, J. P. O striker and N.Y. Gnedin 



The luminosity of accreting BHs in each host galaxy is, by as- 
sumption, proportional to the star formation rate in that galaxy, that 
is, to the ability of the gas in that galaxy to cool efficiently and be 
accreted toward the centre. Due to the early reheating by the X-ray 
background and feedback processes (i.e., galactic winds produced 
by photoionisation and SN explosions) the star formation in the 
smaller mass galaxies is strongly suppressed and so is the accre- 
tion rate on intermediate mass BHs. Moderately larger galaxies are 
able to retain their gas but they cool so inefficiently that few or no 
stars at all are formed. When seed BHs are accreted on these galax- 
ies they find a large reservoir of gas and they can start accreting 
efficiently. 

The X-ray luminosity function, shown in Fig. [9] is dominated 
by the most massive galaxies that are less affected by feedback ef- 
fects and can retain gas more effectively. This means that the most 
rare and massive galaxies at z ~ 15 are expected to be the most lu- 
minous X-ray sources. The volume of our simulation (1 h _1 Mpc) 3 
is large enough to be a representative region of the universe at 
z ~ 15 but it is too small to contain any of these rarer objects. We 
find that the luminosity of the typical X-ray sources at z ~ 15 is 
two orders of magnitude fainter than Chandra dee p field 3cr detec- 
tion li mit for 1.945 x 10 6 seconds integration time I Alexan der et alJ 
l2003t) (limiting flux of 2.5 x 10" 17 (1.4 x 10" 16 ) erg s _1 cm -2 
sr _1 in the 0.5 — 2(2 — 8) keV band). But we expect a few thousands 
of these sources per arcmin 2 . We then use the Press-Schechter for- 
malism to extrapolate the number counts found in our simulation to 
larger volumes. The extrapolated high luminosity tail of the num- 
ber counts is shown in Fig.|5|with a dotted line. The plot shows that 
it is plausible to find a few sources per arcmin 2 above the detec- 
tion limit (shown by the vertical line) of the Chandra X-ray deep 
field (unless the high tail of the luminosity function has a sharper 
cutoff than in the semianalytic predictions). The triangles show the 
number counts of point sources of the 2 Ms Chandra deep field 
in the 0.5-2 keV band and the closed polygon shows the expected 
slope of the luminosity function based o n the statistics of pix el fluc- 
tuation from unresolved point sources iMivaii & Griffiths! £2002). 
The solid line that fits the number counts is a model prediction 
for the expected counts from AGNs (e. g. . IComastri et aHll995l) . 
We see that the observed steepening of the number counts slope at 
fluxes Fx ,$ 10~ 16 is consistent with the predictions of our simu- 
lation but it is also consistent with the predicted X-ray contribution 
by starburst galaxies (e. e. JPtak et alj200lllRanalli et all2003l) . In 
principle it is possible to separate the contribution of these two pop- 
ulations since the high-z sources, contrary to starburst galaxies, do 
not have optical counterparts but may be detectable in the infrared if 
they are not obscured locally by dust. The planed spatial resolution 
of Constellation-X and XEUS of about 5 arcsec is not sufficient to 
avoid source confusion (Mushotzky, private communication). The 
cross correlation of the unresolved X-ray background fluctuation 
with the optical number counts may be used to separate the contri- 
bution of the high redshift population from the starburst population. 

This prediction has exciting consequences for the direct ob- 
servability of X-ray sources at re dshifts as high as z ~ 15. In- 
terestingly, Koe kemoer et alJ 120031) describe a possible new class 
of X-ray sources that have robust detections in ultra-deep Chandra 
data, yet have no optical counterpart in deep multi-band GOODS 
Hubble Space Telescope (HST) ACS images. Their ratios of X-ray 
to optical fluxes are at least an order of magnitude above those gen- 
erally found for other AGN, even those that are harboured by red- 
dened hosts. The authors conclude that if these sources lie above 
redshifts 6-7, such that even their Lyman-a emission is redshifted 
out of the bandpass of the ACS z(850) filter, then their optical and 
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Figure 10. Top panel: the mean redshifted 21 cm signal as a function of 
frequency and redshift for the simulation run M-BH (thin lines) and our 
fiducial semi-analytical model (thick lines). The two hatched zones show 
the locations of TV and FM radio bands within which is likely that no ob- 
servations will be possible due to man-made interference. Bottom panel: 
the fluctuations in the redshifted 21 cm signal in a beam of 10' and 1MHz 
bandwidth. The solid black band is the sensitivity limit of a filled dish aper- 
ture telescope for 1000 hours of integration. 



X-ray fluxes can be accounted for. They find seven of these sources 
in the field, in agreement with our estimates for the M-BH run. In 
our model, if the source emission extends to the rest frame visi- 
ble/UV bands (i.e., if they are not heavily reddened by dust), they 
should be detectable in the infrared above 2-3 micron. Otherwise 
they should be visible only in X-ray bands and in the FIR at about 
100 Atm- 

The postulated existence of this mini-quasar population in- 
creases the probability of finding a few radio loud sources at 
z ~ 17, that could be successfully used to study the redshifted 
21 cm absorption lines produced by the hyperfine transition of neu- 
tral hydrogen in the IGM before reionisation JCarilli et alj|2002l) . 
Next we discuss further the radio signal at 21 (z + 1) cm. 



4 REDSHIFTED 21 CM LINE IN EMISSION AND 
ABSORPTION 

One possible observational test for our model is the de- 
tection of redshifted 21 cm radiation from high redshift 
iMadau. Meiksin. & Reeslll997LlTozzi et alj|2000h . Fig. flOl shows 
the predicted mean and rms signal for our representative simula- 
tion and the best-fit fiducial model. The only instrument capable 
of detecting the redshifted 21 cm signal is the projected Low Fre- 
quency Array (LOFAR). The solid black band shows the expected 
LOFAR sensitivity for observations with the core (2 km in diame- 
ter) for 1000 hour integration time. 

In calculating the expected 21 cm signal, we included three 
physical effects that influence the level population for the hyper- 
fine transition in hydrogen atom: pumping by Lyman-a photons 
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and collisions with electrons and neutral atoms <Tozzi et all200Ch . 
The former is the dominant effect at lower redshifts (z < 20), while 
the latter is important at higher redshifts. We then simulated a num- 
ber of lines of sight through the computational box and averaged 
over them to properly include the effect of velocity focusing: be- 
cause most of the emission comes from the high density regions, 
which are collapsing, velocity focusing increases the effect by a 
small but not negligible factor. 

Fluctuations in the expected 21 cm signal, however, come 
from large spatial scales (larger than about 10 Mpc). The density 
fluctuations in this regime are linear, and can be computed analyti- 
cally (Gnedin & Shaver 2004). 

As one can see from Fig. llQI the mean signal is well within the 
sensitivity of existing or proposed low frequency radio telescopes 
(such as Arecibo, LOFAR, or SKA). However, the main challenge 
to observing the mean signal is not the sensitivity, but rather the 
foreground contamination from the Galaxy (cf , IOh & Mackl2003l) . 
The smoothly variable emission signal will most likely be unob- 
servable, but in our model a remarkable opportunity arises: a sharp 
absorption feature at z ~ 25 falls just outside the broadcast TV 
band (that starts at 54 MHz, which corresponds to z — 25). This 
feature will appear as an about 100 mK "absorption line" in the 
spectrum of the Galactic foreground, and may be observable with 
future radio telescopes. The specific location of this feature is, of 
course, model dependent - a slight variation of cosmological pa- 
rameters well within WMAP errors can either move it to higher 
redshifts, or hide it entirely behind the broadcast TV band. Thus in 
the later case, observing it will probably be impossible, but if in- 
deed it falls beyond the broadcast TV band, it will clearly indicate 
the "first light" in the Universe - or, at least, first Lyman-a photons. 

The observational situation with fluctuations is, in some sense, 
reversed with respect to the mean signal. While it will be signifi- 
cantly easier to separate fluctuations in redshifted 21 cm emission 
from the Galactic foreground, the early emission of Lyman-a pho- 
tons relative to the growth of structure makes the fluctuation signal 
harde r to observe than i n more conventional reionisation models 
iGnedin & Shavei]|2003t) . In addition, at the predicted level of the 
cosmological signal, systematic errors such as beam leakage also 
become important. It, therefore, appears that unless the absorption 
feature at 54 MHz can be detected, existing and projected radio 
telescopes (including LOFAR and SKA) will not be able to detect 
cosmological fluctuations in the redshifted 21 cm signal in emis- 
sion (but there still remain the possibility to detect the 21 cm forest 
in absorption against high redshift radio sources). 



5 SECONDARY CMB ANISOTROPIES 

Another possible observational signature of the early episode 
of X-ray emis s ion is secondary CMB anisotropies. We used 
IGnedin & Jaffd (2001) method to compute the spectrum of sec- 
ondary anisotropies on small angular scales (arcmin 2 or I > 10 4 ), 
where the anisotropies from the first episode of structure formation 
dominate other contributions. Maps of the temperature anisotropies 
for the X-ray preionisation simulation (run M-BH, left panel) and 
Population III stellar reionisation simulation (run 128L2VM in pa- 
per I, right panel) are shown in Fig. 1121 

Fig. 1121 shows the angular spectrum for runs M-B H and 
128L2 VM together with the typical spectrum from lGnedin & Jaffel 
1200 it) . A remarkable feature can be observed: while in a model 
without early episode of X-ray emission Cil(l + 1) falls off 
at small angular scales approximately as l~ 3 ^ A , in our model 
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Figure 11. Secondary temperature anisotropies (in fiK) in a 0.55 X 0.55 
arcmin 2 patch of the sky for the run M-BH. The map shows strongly non- 
linear features produced by gas in minihaloes. 

the anisotropies instead grow roughly as The temperature 
anisotropies are of the order of AT ~ 16/iK at ~ 1 arcsec scales. 
This opens up a unique opportunity to test our models against those 
without early episode of X-ray emission, although, observations on 
these scales are in the future and may be extremely difficult to sepa- 
rate the signal from the foreground sources contamination (but see, 
iFomalont et alll993tlchurch et alll997h . 

The non-linear Ostriker-Vishniac effect (O striker & Vi shniac 
1986) describes anisotropies generated in a universe in which the 
ionisation fraction is homogeneous in space. In our simulations 
the ionisation fraction is not homogeneous in space and the total 
power spectrum of the anisotropies includes a contribution from 
the "patchy reionisation". In the X-ray preionisation model the non- 
linear OV component of the anisotropy power spectrum (calculated 
assuming the mean ionisation fraction) is actually slightly larger 
than the total signal, which means that the "patchy reionisation" 
component is mildly anti-correlated with the OV component. In 
other words, the higher density regions are less ionised than the 
low density regions - this is, obviously, a signature of X-ray ionisa- 
tion. Such a separation is purely artificial and unphysical, but illus- 
trates the main physical difference with the UV-ionisation models. 
Indeed. IGnedin & Jaffel l200ll) have shown that in models in which 
reionisation is produced by UV from stellar sources the OV and 
"patchy reionisation" signals are instead strongly correlated. 

The temperature anisotropies on these scales are dominated 
by non-linear structures and an analytical derivation is not possible. 
But it is clear from Fig.[2|that the early start of IGM ionisation in 
the X-ray models produces additional power on arcsec scales that 
should be roughly proportional to the mean fractional ionisation 
of the IGM during the extended period of X-ray pre-ionisation. In 
principle, the detection of anisotropies on these scales can be used 
to obtain precise measurements of the ionisation history of the IGM 
before the epoch of overlap when the reionisation is competed. 



6 DISCUSSION AND SUMMARY 

The present paper is the third in a series (see paper I and paper Ha) 
devoted to the study of physically plausible reionisation scenarios. 
The models have optical depths to Thompson scattering r c ~ 0.17, 
as measured by WMAP ( K ogut et all2003l) . and are consistent with 
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Figure 12. The power spectrum of the CMB anisotropics for the run M- 
BH (solid thick line). For comparison, the dashed and dotted lines show 
the power spectrum for two models without early X-ray sources: the dashed 
lines is the simulation 128L2noSN with box size of 2 h" 1 Mpc (from Ri- 
cotti & Ostriker 2003) and the dotted line is a simulation with box size of 4 
h Mpc (from Gnedin & Jaffe 2000). Two thin straight lines mark +0.75 
and -0.75 slopes respectively. Two vertical arrows show the fundamental 
and the Nyquist frequencies of the simulated sky for the run M-BH. 



observations of the IGM optical depth to Lyman -q an d Lyman- 
pho tons toward high redshift quasars iFan et al]|2003t Isongaila 
2004). We have studied models of reionisation by stars (paper I) 
or mini-quasars (paper Ha) considering realistic physical scenar- 
ios and observational constraints. Our main conclusion is that it 
is very difficult for any of these models to produce optical depths 
in excess of 0.17. If this large optical depth is produced by ultra- 
violet radiation from stellar sources (i.e., Population III stars with 
a top-heavy IMF), we find that zero-metallicity stars must be the 
dominant mode of star formation up to redshift z ~ 10. If this 
scenario is correct, we need to understand the reasons for the very 
inefficient mixing of metal enriched gas from SNe with the gas in 
which star formation takes place. At the moment this inefficient 
mixing is not reproduced by numerical simulations. A feasible al- 
ternative scenario requires that most Population III stars must im- 
plode into black holes without exploding as SNe or have sublumi- 
nous SN explosions with a large amount of metal fall-back onto the 
compact remnant in order to reduce their metal yields. A scenario 
in which a large fraction of Population III stars end their lives as 
pair-instability supernovae is not compatible with the large optical 
depth to Thompson scattering measured by WMAP 

In this paper we use cosmological simulations to study a reion- 
isation scenario in which standard reionisation by Population II 
stars is preceded by partial ionisation and reheating at early times 
by an X-ray background (i.e., "X-ray preionisation" models stud- 
ied in paper Ha using semianalytic simulations). In these models the 
ionisation rate from the secondary radiation produced by accretion 
on compact remnants from the first stars dominates over the pri- 
mary ultraviolet radiation emitted by the stars during their lifetime. 
The most appealing aspect of these models is their insensitivity to 
the duration of the epoch of Population III stars domination. In fact, 
the relative importance of Population III stars with respect to Popu- 
lation II stars in practice has to be treated as a free parameter due to 
the large uncertainties in modelling the complex physical processes 
that regulate metal production and mixing. 

In order to reproduce WMAP results we find that a frac- 
tion of about 10~ 4 of all the baryons in the Universe needs to 



be accreted by compact objects before redshift 10. This fraction 
is comparable in mass to the observed mass density of supermas- 
sive BHs in the galactic nuclei today. This does not pose a strong 
constraint on the models because a sizeable fraction of these early 
intermediate-mass BHs is expected to be removed from the galax- 
ies during the last phases of the merger or to reside in the inter- 
stellar medium of galaxies, being hardly detectable lMadau & Reesl 
(2001). Perhaps a small fraction of these intermediate-mass BHs 
may be accreting and contribute to the obser ved population of ULX 
lAgol & Karmonkowskj200llKrolikl2'0bl) . 

The stronger constraint on the model is posed by the observed 
soft X-ray background. Assuming that the early mini-quasars popu- 
lation produces an optical depth to Thompson scattering r c ~ 0.17, 
their contribution to the background in the 5-50 keV bands is 5- 
20%. Future X-ray missions may be able to detect these sources 
if they exist. The predicted fluxes of the most rare objects is at 
the limit of the detection of the Chandra deep fields and about 
1000 objects with flux 10~ 18 erg cm~ 2 s _1 should be present per 
5.3 arcmin 2 . Note that these calculations assume a spectral energy 
distribution of the sources that has a non-negligible high energy 
power-law component (cf., Fig. 4 in paper Ha), in agreement with 
observations of QSOs and ULX in the local universe. If the source 
spectra are dominated by a multicolour disk thermal component 
their contribution to both the X-ray background in the 2-50 keV 
bands and the faint source counts would be negligible. 

The redshifted X-ray background also has interesting conse- 
quences for the reionisation history of He and the thermal history 
of the IGM at redshift z ~ 3, independently of the assumed spec- 
tra of the sources. In this paper we confirm the results of paper Ha 
in which we found that He II is almost fully reionised for the first 
time at redshift z ~ 17 and afterwards slowly recombines before 
experiencing a second reionisation at redshift z ~ 3 produced by 
the redshifted X-ray background. The heating rate from the back- 
ground radiation keeps the temperature of the IGM at about 10,000 
K, in rough agreement with observations of the line widths of the 
Lyman-a forest at z ~ 3 — 4. We also emphasise that the red- 
shifted X-ray background is important in producing a fully ionised 
atomic hydrogen in the low density intergalactic medium before 
stellar reionisation at redshift z ~ 6 — 7. As a result stellar reioni- 
sation is characterised by an almost instantaneous "overlap phase" 
of H II regions. 

The patchy topology of reionisation produced by stellar 
sources contrast with the spatially homogeneous partial ionisa- 
tion by X-rays. This produces distinctive signatures on temper- 
ature/polarisation CMB anisotropies and on the redshifted 21cm 
emission/absorption from neutral hydrogen at high redshift. 

(i) The power spectrum of the EE polarisation is sensitive to the 
visibility function g(z), defined in equation The Plank satellite 
should be able to distinguish between the visibility function pro- 
duced by an early X-ray partial ionisation (cf., Fig. [4] and Fig. 8 
in paper Ha) or the one expected for reionisation by stellar sources 
(cf., Fig. 4 in paper I). 

(ii) On small angular scales (< larcmin or I > 10 4 ) the sec- 
ondary anisotropies produced by non-linear structures during the 
early reionisation epochs dominate over other contributions, offer- 
ing a unique opportunity to study the first episode of structure for- 
mation. A remarkable feature can be observed: while in a model 
without early episode of X-ray emission the power spectrum falls 
off at small angular scales approximately as i -3 / 4 , in X-ray preion- 
isation models the power instead grows roug hly as l 3/i . The tem- 
perature anisotropies are of the order of AT ~ 16/jK at ~ 1 
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arcsec scales. But it is extremely difficult to detect anisotropies 
on t hese scales, partly because of fore ground sources contamina- 
tion iFomalont et alll993tlChurch et alll997l) . In models in which 
reionisation is produced by UV from stellar sources the secondary 
anisotropies produced by the nonlinear Ostriker-Vishniac effect 
(i.e., assuming uniform ionisation fraction) have less power than 
in the full calculation that includes the "patchy reionisation" signal 
(the signals are correlated). In X-ray preionisation models the non- 
linear Ostriker-Vishniac effect and "patchy reionisation" signals are 
instead slightly anti-correlated. 

(iii) The redshifted 21cm emission/absorption from neutral hy- 
drogen is another powerful probe of the ionisation state of the IGM 
at high redshift. In X-ray preionisation models a sharp absorption 
feature at z ~ 25 — 30 falls just outside the broadcast TV bands. 
This feature will appear as a roughly 100 mK "absorption line" 
in the spectrum of the Galactic foreground, and may be observable 
with LOFAR. The partial ionisation of the IGM in the X-ray preion- 
isation scenarios offers better opportunities to observe the fluctua- 
tions of the redshifted 21cm line in emission at redshifts lower than 
10, where the signal is easier to measure. 

Finally indirect signatures of X-ray preionisation are related 
to the discovery of massive BHs in the nuclei of dwarf galaxies 
and/or the identification of intermediate mass BHs in the ISM of 
galaxies (e.g., ULX). If Population III stars are instead responsible 
for the large optical depth measured by WMAP future observations 
with IWST and ground based near-infrared Lyman-a surveys us- 
ing gravitational lenses may soon be able to observe these objects 
IPello et alj 2004: Ricotti et al. 2004) and probe the ionisation state 
of the IGM. 
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